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Abstract 

A two-dimensional Frenkel-Kontorova model under a steady exter- 
nal force is used to study the nonlinear sliding friction between flat 
macroscopic surfaces with a lubricant layer in between. The nonequi- 
librium properties of the model are simulated by a Brownian molecular 
dynamics and results are obtained as a function of temperature and a 
microscopic friction parameter 7]. 
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1 Introduction 



The behavior of an adsorbed layer under an external force in the plane is 
a fundamental problem for understanding the mechanisms responsible for 
friction and lubrication between two flat macroscopic surfaces |l|, 0. Because 
of the lubricant, the friction between the surfaces is essentially determined by 
the force required to shear the adsorbed layer. For small external force, this is 
just the familiar study of the collective diffusion of the overlayer. However, for 
strong external fields, even a pinned layer can acquire significant mobility at 
low temperatures due to the lowering of the diffusion barrier by the external 
field. Moreover, the steady state effective temperature of the overlayer can 
be significantly higher than that of the substrate. Thus we can actually be 
dealing with the flow of a liquid layer rather than simple diffusion in the 
solid phase. Recently, there have been a great deal of interest in the study of 
this problem in regards to boundary lubrication. In particular, Persson 
has made a series of studies based on a model of adatoms interacting with 
Lennard- Jones pair interaction. It is found that the response of the overlayer 
can develop hysteresis. This means that the threshold field Fa required to 
initiate the motion of the adlayer is significantly larger than the threshold 
force Fh below which an initially sliding layer will be pinned again. In fact, 
based on general hydrodynamical considerations, it is argued that the ratio 
of Fa to Fb should be 2 for all systems with sizable frictional coupling to 
the substrate. This hypothesis was tested for the Lennard Jones system by 
direct numerical simulation and seem to be satisfied |^. The universality of 
this relation between sliding friction and static friction still remains an open 
question. Even the existence of a hysteresis effect is not necessarily true for 
other models of adsorbed layers. 

We introduce a simple model for sliding friction based on a two dimen- 
sional version of the Frenkel-Kontorova model Q , where particles are coupled 
to each other by elastic springs and to a periodic potential. The equilibrium 
properties of this model are well understood both in one and two dimen- 
sions, largely in relation to the study of commensurate and incommensurate 
transitions 0] . It has also been a standard model in the study of mobility of 
charge density waves . The infinite range of the harmonic interaction in the 
model leads to a disordered (fluid) state very different from the conventional 
one where atoms interact via short range forces. Even in the fluid phase, 
the atoms cannot diffuse around and will maintain their original neighbours. 
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Thus it is of particular interest to test whether the universahty of hysteresis 
and the relation between static and kinetic friction still applies to this simple 
model. 

2 Model and Simulation 

The simplest version of the Frenkel-Kontorova model in two dimensions is to 
limit the displacements of the individual atoms in one direction only. Thus 
the Hamiltonian is of the form 

— t/o cos(27rxjj/a)} (1) 

which can be taken to represent a system of coupled linear harmonic chains 
of particles, each of mass m moving in a periodic potential of period a and 
amplitude Uo- The second term represents the elastic interaction within a 
chain and the third term represents interaction between adatoms on different 
chains. The index i,j labels the particle with integer i {i = 1, ) and the 
chain with integer j {j = 1, ...,Ny), so that the system consists of A^^ x A^^ 
adsorbates at positions Xij. We have considered = Ny = N and = 

Ky = K. 

Most of our results are obtained assuming periodic boundary conditions 
which means that the N + l,j particle is a periodic image of particle 1, j in 
each chain, so that Xjv+i.j = xij + A^^ c^, where A^^ is the number of local 
minima of the periodic substrate in the x-direction. In the other direction, 
periodic boundary condition means that Xi^N+i = Xi^i. The average distance 
between particles, b, in each chain is determined through the boundary con- 
dition and is given hj b = Ns a/N and the overlayer coverage can be defined 
as 6* = N/Ng. With this setup, b and 9 cannot be changed independently 
and the ratio b/a is forced to be rational. We have also performed few cal- 
culations with free boundary conditions in the x-directions. In this case, the 
ratio b/a can assume arbitrary nonrational values. The system now can re- 
main commensurate for a given 6 only up to a critical value 6c of the misfit 
parameter [Q] 6 = {b — a)/a. 
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To study sliding friction in the model we assume that the substrate acts 
as a heat bath at an equilibrium temperature T and that an external force 
F acts on each of the adsorbates. The equation of motion for the particle 
coordinates are described by Brownian dynamics 

dV dU ^ ^ 
mxij + mriXij = \- F + fij (2) 

kJ ilj 'ij kJ ilj 'Ij 

where U is the periodic potential, V is the adsorbate-adsorbate harmonic in- 
teraction potential and fij is a random force that is related to the microscopic 
friction parameter t] by the fluctuation dissipation theorem 

< fij{t)fi'j'{t') >= 2r]mkBT {t - t') (3) 

We have studied the sliding friction of the adsorbate by simulating the 
above equations using standard methods of Brownian molecular dynamics [^]. 
We use units in which a = 1, m = 1 and Uo = 1- Typically, the time variable 
was discretized with time step 5t = 0.02r to 0.06r where r = {ma^ /UoY^'^ 
and 4 x 10^ time steps were used in each calculation allowing 10^ time steps 
for equilibration. The elastic constant was set to = 10 and calculations 
were performed as function of temperature, external force and coverage for a 
system containing N x N adsorbate particles where = 10. The system was 
allowed to evolve into a steady state such that the time average of physical 
quantities approached a constant. 



3 Results and Discussion 

When the external F = 0, the drift velocity Vd of the adsorbate center of 
mass is zero independent of the temperature with each particle performing 
a Brownian motion. At T = OK, the velocity remains zero for any applied 
force less than a critical value Fg where the effective potential minima disap- 
pears, Fo = 27r[/o, the same result as for the one-dimensional case 0]. From 
now on, we will normalize the external force by this critical value. If a small 
force is applied at higher enough temperature, in addition to the irregular 
motion of the particles, one expects an overall drift in the direction of the 
force with a speed va proportional to the external force, Vd = F/mfj in the 
linear regime, where f] is an effective sliding friction. Fig. 1 shows the results 
of the inverse of the normalized linear sliding friction, rj/fj, calculated as a 
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function of temperature for a small value of the external force, F/ Fo = 0.035 
and b/a = 2. At low temperatures the adlayer is in a pinned state at small 
forces, and rj/f] is essentially zero. As the temperature is increased 1/f/ re- 
mains zero until T = T^, beyond which T]/f] increases significantly and should 
approach 1 at higher enough T. The value of correspond to the commensu- 
rate solid-fluid transition where the overlayer at high temperatures become 
disordered and is effectively depinned from the substrate. Figs. 2 and 3 
shows snapshot pictures of the adsorbate in the commensurate solid and the 
fluid state respectively. From Fig. 1 we estimate that Tc ~ 2.3. We have also 
performed similar calculations with different boundary conditions. In Fig. 4 
we show the results for the structure factor at low and high temperatures for 
b/a = 1.9 with free boundary condition. At low T the structure factor has a 
sharp peak at g = vr and the overlayer is commensurate with the substrate 
with lattice vector 2a. At high T the structure factor is much broader corre- 
sponding to a fluid state. For b/a = 1.9 it peaks at an incommensurate value 
which is temperature dependent and limited by g = 27r/6. The corresponding 
result for b/a = 2 and p.b.c is similar except that the structure factor at high 
temperature remains peaked at g = tt. 

In the nonlinear regime, the relation between Vd and F is expected to 
depend strongly the initial phase [||. If the adsorbate, when F = 0, is a 
fluid, Vd will be nonzero for arbitrarily small external force. This result is 
also true for the present model. In the fluid phase above Tc, the overlayer 
is essentially depinned and any small force leads to nonzero Vd- Moreover, 
the relation Vd = f{F) shows no hysteresis. For T < T^. however, where 
a commensurate (solid phase) prevails when F = 0, a critical force Fa is 
necessary to depin and initiate sliding. It is expected [0] that in this case 
the nonlinear sliding friction exhibit hysteresis as a function of F, i.e., the 
relation between velocity and external force depends on whether F increases 
from zero or decreases from a high value and this in turn implies a stick and 
slip motion. 

In Fig. 5, we show the results for Vd as a function of F for r] = 0.6, 
T = 0.5. Note that in this case the temperature is considerably less than 
the activation energy barrier {2Uo) for particle diffusion. We found that this 
is an important requirement for hysteresis to appear. Initially, the adlayer 
is in a pinned state. Then the simulation was started with various external 
force strengths. An hysteresis loop is clearly seen, corresponding to a stick 
and slip motion. The effective temperature T* of the adlayer is shown Fig. 
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6. It is seen that the initial shding phase when the apphed force is increased 
beyond the static threshold Fa has almost the same temperature as the sub- 
strate and the substrate potential provides no additional resistance. This is 
a sliding solid phase. It corresponds to the floating or incommensurate solid 
phase in the equilibrium situation. It is important to note that this phase 
is dynamically generated by the external field. With our present choice of 
parameters, the misfit is zero and the floating solid phase does not exist in 
equilibrium. Thus the external force induces a commensurate solid to sliding 
solid transition aX F = Fa- As the force is decreased through F^, there is 
a jump in T* accompanied by a simultaneous decrease in Vd- This corre- 
sponds to an overheating of overlayer to an effective temperature above 
(compare with Fig. 1), dynamically "melting" the sliding solid phase. This 
can be regarded as sliding solid to fluid transition. Finally at F = Ff,, the 
fluid phase condenses into a commensurate solid phase and the tempera- 
ture of the adlayer and the substrate become the same again. The fact that 
there is a condensation from the "melted" phase is confirmed by noting that 
the temperature of the adlayer at this condensation is almost identical to 
the temperature Tc, the temperature where the pinned solid to fluid phase 
transition occurs. Thus there is already a rich variety of dynamical phase 
transitions manifested in this simple model. We note that Persson p| has 
observed very similar phase transitions for Lennard- Jones systems, showing 
that these features are universal in the non-linear response of overlayers and 
not just for specific models. In addition the present results show that a con- 
ventional "liquid" phase is not required and a high temperature disordered 
fluid phase with vanishing shear resistance is sufficient to show some of these 
effects. 

We now consider the behavior of the nonlinear sliding friction at higher 
temperatures but still less than Tc. We flnd that at T ^ f/o, hysteresis 
disappears as indicated in Fig. 7 but Fa can still be deflned. This behavior 
is not speciflc to the present model and has also been found in the molecular 
dynamics simulation of the Lennard- Jones system 0. This temperature is 
comparable to the energy barrier for center of the mass diffusion and for 
a given force F, within the range where hysteresis is observed at low T, 
the system is able to relax in a short enough time to a stable single velocity, 
independent of the history of force variation. This suggests that the hysteresis 
observed at low T is due to metastability. 

We have also studied the effects of varying the microscopic friction pa- 
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rameter 77 and found that the ratio Fh/Fa strongly depends on rj and converge 
to Fh/ Fa = 1 for sufficiently large 77 as shown in Fig. 8. This should be con- 
trasted to the results obtained by Persson [0 for the Lennard- Jones system 
where F^/Fa is found to saturate for large r] at a universal value ~ 0.6. We 
suggest that this is due to the fact that the nature of the fluid phase in the 
two models are very differrent. The hydrodynamics arguments used by Pers- 
son based on the existence of a drag force acting on a nucleating domain 
may not be valid in this case. 

4 Conclusions 

There are many questions that remains to be investigated. The most impor- 
tant one is probably to establish which features of the nonlinear transport 
are universal and model independent and which ones are model sensitive. 
Already, we have observed in the Frenkel-Kontorova model that under in- 
creasing values of the microscopic friction rj, the hysteresis loop tends to 
disappear whereas it is stable in the Lennard Jones system. We suspect that 
this is related to the stability of the sliding solid phase as a function of rj 
as well as the dependence of the fluid-pinned solid transition on the nature 
of the fluid state. However, many features observed in the Lennard- Jones 
system also appear in the simple Frenkel-Kontorova model. The advantage 
of this model system is that the equilibrium property of the commensurate- 
incommensurate transition is well understood through a series of studies. 

t Research supported in part by a joint NSF-CNPq grant. 
* Research supported in part by an ONR grant. 
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Figure 1: Inverse of the normalized sliding friction 77 as a function of tem- 
perature T ior rj — 0.5, b/a — 2 
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Figure 2: Snapshot picture of the adsorbate in the commensurate sohd ph; 
for T = 0.5, b/a = 2. 
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Figure 3: Snapshot picture of the adsorbate in the high temperature phase, 
for T = 3, h/a = 2. 
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Figure 4: Structure factor S{q) as a function of reciprocal lattice vector q 
at different temperatures, using free boundary conditions and b/a = 1.9. 
Open diamonds correspond to T = 0.5 and circles to T = 3. The large and 
small arrows indicate reciprocal lattice vectors q = 27r/2a and q = 27r/6, 
respectively. 
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Figure 5: Drif velocity Vd as a function of external force F for rj = 0.6 and 
r = 0.5. 
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Figure 6: Effective temperature T* of the overlayer normalized to the sub- 
strate temperature T as a function of the external force F ior r] — 0.6 and 
r = 0.5. 
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Figure 7: Drif velocity Vd as a function of external force F for T = 1 and 
77 = 0.5. 
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Figure 8: Ratio Fi,/Fa as a function of 77 for T = 0.5 
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